Nonlinear unbalanced Bessel beams in the collapse of Gaussian beams 

arrested by nonlinear losses 
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Collapse of a Gaussian beam in self-focusing Kerr media arrested by nonlinear losses may lead 
to the spontaneous formation of a quasi-stationary nonlinear unbalanced Bessel beam with finite 
energy, which can propagate without significant distortion over tens of diffraction lengths, and 
without peak intensity attenuation while the beam power is drastically diminishing. 



Nonlinear unbalanced Bessel beams (NL-UBB) [3| 
have been recently described as the only possible type of 
light beam that can propagate without distortion and at- 
tenuation in a medium with nonlinear mechanisms of dis- 
sipation of energy. Contrary to spatial solitons, NL-UBB 
transport infinite power, which allows for a continuous 
transversal energy flux from the beam periphery towards 
its nonlinear center, refuelling the power absorbed in this 
region during propagation. [l|] The idea of NL-UBB has 
been extended to Bose-Einstein condensates, where simi- 
lar long-lived structures are supported by the dissipation 
provided by three-body inelastic collisions. Q 

A spatial dynamics showing the spontaneous trans- 
formation of a Bessel beam into a NL-UBB in a Kerr 
medium with nonlinear losses (NLL) has been observed 
in Q ■ More importantly, the transformation of an ultra- 
short pulsed Bessel beam into a pulsed NL-UBB may take 

glace without the excitation of any temporal dynamics, 
I as the spectral broadening and temporal splitting that 
would occur in the Bessel beam in absence of NLL. These 
NL-UBB stabilized by NLL have indeed been observed in 
recent experiments. t4| 

The robustness of the NL-UBB suggests that it may 
act as an attractor of the Kerr dynamics with NLL over 
a broader space of input conditions. In this Letter we 
report on numerical simulations of the propagation of an 
input Gaussian beam showing the spontaneous formation 
of a quasi-stationary light beam with most of relevant fea- 
tures of a NL-UBB. The different regimes of propagation 
of a Gaussian beam in a medium with Kerr nonlinearity 
and NLL were identified in [H . NL-UBBs are seen here to 
emerge in case of strong self-focusing stopped by a small 
amount of NLL. Contrarily to a Bessel beam, the input 
Gaussian beam carries finite power. As a consequence, 
only the nonlinear core of the NL-UBB is formed. Sta- 
tionarity of the nonlinear core is seen to be sustained by 
an energy flux from the beam periphery, as in NL-UBBs. 
However, since the available power is finite and decreas- 
ing, stationarity lasts after a propagation distance, which 
is, nevertheless, much larger than the diffraction length 
associated to the nonlinear core. Our monochromatic 
approach is expected to describe also experiments with 
short pulses insofar NLL is the mechanism stopping col- 
lapse, and can be assumed to quench any significant tem- 



poral dynamics afterwards. [4j 

The nonlinear Schrodinger equation (NLSE), [J, [3] 



2k r n 



K) 



-\A\ 



2K-2 



A, 
(1) 

is assumed to describe the propagation of a light beam 
E = A(r, z) exp(— ia;t -t- ikz) with revolution symme- 
try along the z axis in a medium with Kerr nonlin- 
earity (linear and nonlinear refraction indexes n and 
n2 > 0, respectively) and multi-photon absorption (coef- 
ficient > 0, if = 2,3...). For the incident Gaus- 
sian beam A{r,0) = \/7o exp(— r^/sg) of peak intensity 
lo and width so, we introduce dimensionless amplitude 
A = A/^/Io, radial coordinate p = r/so, and propagation 
distance ^ = z/ksg, to write ^ as 
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(a) On-axis intensity |j4(0, ^)p for input Gaussian 
: 60 and increasing 7 {K = 8). (b) Intensity |j4(p,^)p 



for g = 60 and 7 = lO"'' {K - 



(in logarithmic scale). 
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FIG. 2: Beam power 27r JJ^ dpp|j4(p,^)|^ as a function of 
propagation distance ^ for g = 60 and increasing 7 {K = 8). 



and the initial condition as A{p,Q) = exp(— p^), where 
g = k'^n2losl/n and 7 = P^^^ lQ~^ksl/2 measure the 
strengths of Kerr nonhnearity and of NLL in the input 
Gaussian beam. 

Without NLL, the Gaussian beam collapses if g > 
go ~ 3.79. Figure [1] (a) displays the peak intensity 
/ = |A(0,^)p along ^, numerically evaluated from ([2|) 
with 5 = 60 ^ where the strong collapse is arrested 
by different amounts of NLL 7 {K — 8). Similar trend 
is observed for other values of g, 7 and K such that the 
characteristic lengths ^sf — l/(2y^), ■Co — 1/2, and 
^NLL = (2^-^ - l)/[27(i^ - 1)] of self-focusing, diffrac- 
tion and NLL of the input Gaussian (defined as in [^) 
satisfy ^sf <C ^ ^nll- Under these conditions NLL 
are negligible in the fast self-focusing, becoming impor- 
tant only in the collapse region. Arrest of collapse by 
NLL leads to an axial region of sustained high intensity. 
A rough estimate of the peak intensities / (in units of the 
input intensity Iq) involved in this region is provided by 
the simple formula / « (ff/7)^/^^"^^ (e.g., / = 19.8,6.3 
and 2.9 for 7 — 10^^, 10"'^ and 10~^, respectively) which 
is obtained, remarkably, by equating the strength gl of 
the Kerr nonhnearity (for a field of intensity /) to the 
strength 7/^^^ of NLL. With increasing 7, ripples in 
the peak intensity / disappear, the intensity then remain- 
ing nearly constant for intermediate values of 7 (between 
10^"' and 10"'^ for 7 — 60), or monotonically, though 
slowly decreasing for higher 7. 

The transversal profile after collapse features a central 
spike surrounded by one or more rings, as seen in Fig. 
[Ijb) . The width of the central spike can be characterized 
by 1/ \/2gI, and remains therefore constant upon propa- 
gation as long as / is constant. Fig. [T]^b) corresponds to 
the case of flattest peak intensity / (7 = 10~^), and ev- 
idences that the inner portion of the transversal profile, 
comprising not only the central spike but also one ring, 
remains quasi-stationary in the region of flat peak inten- 
sity, while outer rings are seen to spread due to diffrac- 
tion. The axial length ~ 0.05 of this region is estimated 
to be about 17 times the diffraction length s^/2 ~ 0.003 
associated to the width s ~ 0.077 (1/e^ intensity decay) 
of the central spike. 



Figure ([2]) shows that these properties of stationarity 
are accompanied, paradoxically, by a drastic diminution 
of the beam power due to NLL, which are greatly en- 
hanced in this region of high intensity. Being nonlinear, 
absorption takes place mainly in the inner, hot core of 
the beam, just that remaining stationary. Along the re- 
gion of flat peak intensity (7 — 10"'^), for instance, the 
power is seen to diminish down to about one half its ini- 
tial value. 

Stationarity with non-negligible NLL is the distin- 
guishing property of NL-UBBs. 1] These are solutions to 
the NLSE @ of the form A{p, = a{p) exp[i(f>{p) - iSS,], 
with (5 > 0, whose stationary real amplitude a{p) > 
and phase (f>{p) must satisfy, from 121), the ordinary dif- 
ferential equations 
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Eq. Q is a refilling condition, establishing that the 
nonlinear power loss (per unit propagation length) in 
any disk of radius p must equal to the power entering 
into it across its boundary (per unit length) for sta- 
tionarity to be possible. Given 5, 7 and K , NL-UBBs 
solutions to ([3]) and ([U exist with any peak intensity 
/ = a^(0) E^nd axial wave number shift 5 satisfying 
5 ^ min{0, 277i<-7/^~^ — g/}, as illustrated in Fig. [3ja), 
where rjK = 1.67,0.27,0.19,0.16,0.14,0.12,0.11,... for 
K = 2,2), . . .. Far from its nonlinear core, the NL-UBB 
approaches zero as p — > 00 in the form of an unbal- 
anced Bessel beam, a generalization of the well-known 
Bessel beam with unequal amplitudes of its inward and 
outward Hankel beam components, and that carries, as 
Bessel beams, infinite power, [ij Figure[3Ub) shows radial 
intensity profiles and inward radial energy fluxes of NL- 
UBBs with fixed intensity / and decreasing 5, evaluated 
numerically from Eqs. ([3]) and ([4]). At fixed g, 7 and 
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FIG. 3: (a) For g = 60, 7 = 10"^ and K = region of 
existence of NL-UBBs in their parameter space {1,5). (b) 
Radial intensity profiles a?{p) and inward radial energy fluxes 
-2-Kpa^{p)d4){p)/dp for NL-UBBs with fixed I = 6.575 and 
decreasing 5 down to 0"*", as indicated by the arrow in (a)]. 
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FIG. 4; For g = 60, 7 = 10"\ K ^ 8 and input Gaussian 
(dotted curve), radial intensity profiles of the propagated field 
at increasing distances beyond collapse (solid curves), and ra- 
dial intensity profiles of NL-UBBs with same peak intensities 
/ = 3.15, 2.82, 2.23 as those of the propagated field at each 
distance, and 6 = (dashed curves). 

NL-UBB profiles depend on the two parameters / and 5, 
but as 6 diminishes down to (5 = 0+ with fixed / [arrow 
in Fig. [3l^a)], the radial intensity profile of the NL-UBB 
is seen to reach a limiting profile that depends only on 
the choice of / [solid thick curve in Fig. (Sib)]. 

Coming back to the collapse of the Gaussian beam, 
the slow evolution along ^ after collapse is interpretable 
as an adiabatic sequence of NL-UBBs with slowly varying 
peak intensity / and (5 = 0"*". Figure 2] shows, for g = 60, 
7 = 10^^ [lowest curve in Fig. [IJa)], that the inner part 
of the radial intensity profile at each propagation dis- 
tance, comprising the central maximum and about one 
ring, fits well to that of the NL-UBBs with (5 = 0+ and 
with / equal to that of the propagated field. A truly sta- 
tionary NL-UBB cannot be completely formed because of 
the finite power available, which forces the slow evolution 
in the parameter space of NL-UBBs, with peak intensities 
about / = (5/7)^/(^-2) ^ 2.9. The fact that NL-UBBs 
with (5 = 0"*" are always formed can be understood by 
thinking on the collapsing beam as a bundle of rays at 
different angles down to zero, or equivalently, different 
Bessel beams with decreasing 6 (or cone angle). Each 
Bessel beam with 6 Q forms indeed its own NL-UBB 
with same (5 ^ in the nonlinear medium, as demon- 
strated in but only those with 5 — > 0, whose radial 
profile is stationary against a change of S, emerge. 

Figure EJa) is similar to Fig. 2] but for 7 = lO^'^ 
with flatter peak intensity [see also Figs. [Ua) and 
(b)]. The mechanism of energy refilling that sustains 
quasi-stationarity is evidenced in Fig. Eljb). Writing 
^(PjC) = c'CPi 6xp[*'/'(Pi 0]; where a > and (/? are 
the real amplitude and phase, the NLSE ([2|) yields 

—2tt / dppa'^ = -2'Kpa'^dpip - 2727r / dppa^^ , (5) 
Jo Jo 

or dPp/d^ — Fp — Np for short, which generalizes ^ to 
a non-stationary field. Relation (O expresses that the 
variation of the power Pp in a disk of radius p can be 



due to an inward energy flux Fp across the disk bound- 
ary, and to the nonlinear power loss Np in the disk. For 
NL-UBBs, dPp/d^ = 0, and ^ reduces to the refiUing 
condition i.e., Fp — Np. In particular, the total non- 
linear power loss TVqo is compensated in the NL-UBB by a 
constant energy flux Foe coming from infinity [horizontal 
asymptota labelled as NL-UBB in Fig. [H^b)], this being 
made feasible by the infinite power in the NL-UBB. In 
the collapsed Gaussian beam, instead, the total nonlinear 
power loss No^ (gray horizontal asymptotas) is not com- 
pensated by any inward flux from infinity, since F^a = 
(black curves approaching zero) , and the beam power di- 
minishes according to dPoo/d^ = N^o (as in Fig. 
Nevertheless, the collapsed Gaussian beam is seen in Fig. 
[Sfb) to develop an inward energy flux in its inner portion 
such that the refilling condition Fp = Np characteristic of 
the NL-UBBs, and hence dPp/d£_ — 0, are approximately 
satisfied in the same radial region where the radial profile 
propagates undistorted. 

We can then conclude that NL-UBBs tend to be 
formed spontaneously in self-focusing media with NLL, 
not only from a beam with infinite power as a Bessel 
beam, but also from a Gaussian beam that experiences 
strong self-focusing. Due to the finite power available, 
the attracting NL-UBB is never reached completely, and 
is hence hard to be characterized. Nevertheless, the 
numerical simulations indicate that the system seeks 
the NL-UBB with S — O'^ and a peak intensity / ~ 
((7/7) ^/(-^-^^ balancing the strengths of Kerr nonlinearity 
and NLL. 
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FIG. 5; For g = 60, 7 = 10"^, K ^ 8, (a) radial intensity 
profiles of input Gaussian, propagated field and NL-UBB with 
same peak intensity / = 6.61 as the propagated field, and 
(5 = 0"'". (b) Inward radial energy flux Fp (black curves) and 
NLL Np (gray curves) at increasing propagation distances ^. 
For the NL-UBB, inward flux and NLL are equal at any p. 
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